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Abstract 

We investigate the shear viscosity of neutron star matter in the presence of an antikaon con- 
densate. The electron and muon number densities are reduced due to the appearance of a K~ 
condensate in neutron star matter, whereas the proton number density increases. Consequently 
the shear viscosity due to scatterings of electrons and muons with themselves and protons is lowered 
compared to the case without the condensate. On the other hand, the contribution of proton-proton 
collisions to the proton shear viscosity through electromagnetic and strong interactions, becomes 
important and comparable to the neutron shear viscosity. 
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I. INTRODUCTION 



Shear viscosity plays important roles in neutron star physics. It might damp the r-mode 
instability below the temperature ~ 10 8 K [1] . The knowledge of shear viscosity is essential 
in understanding pulsar glitches and free precession of neutron stars [2|. The calculation of 
the neutron shear viscosity {j] n ) for nonsuperfluid matter using free-space nucleon-nucleon 
scattering data was first done by Flowers and Itoh 0,0]. Cutler and Lindblom 5] fitted the 
results of Flowers and Itoh |4( for the study of viscous damping of oscillations in neutron 
stars. Recently the neutron shear viscosity of pure neutron matter has been investigated in 
a self-consistent way [6j. 

It was noted that electrons, the lightest charged particles and neutrons, the most abun- 
dant particles in neutron star matter contribute significantly to the total shear viscosity. 
Flowers and Itoh found that the neutron viscosity was larger than the combined viscosity of 
electrons and muons (rj eil ) in non-superfluid matter jj] . Further Cutler and Lindblom argued 
ffat the electron viscosity was larger than the neutron viscosity in a superfluid neutron star 

1 as Yakovlev and his collaborator showed 
7| • In the latter calculation, the effects of 



Later Andersson and his collaborators as we 
Vefi > Vn in the presence of proton superfluidity 



the exchange of transverse plasmons in the collisions of charged particles were included and 
it lowered the rj etl compared with the case when only longitudinal plasmons were considered 

Q. 

So far, all of those calculations of shear viscosity were done in neutron star matter com- 
posed of neutrons, protons, electrons and muons. However, exotic forms of matter such 
as hyperon or antikaon condensed matter might appear in the interior of neutron stars. 
Negatively charged hyperons or a K~ condensate could affect the electron shear viscosity 
appreciably. 

Here we focus on the role of K~ meson condensates on the shear viscosity. No calculation 
of shear viscosity involving antikaon condensation has been carried out so far. This motivates 
us to investigate the shear viscosity in the presence of an antikaon condensate. The K~ 
condensate appears at 2-3 times the normal nuclear matter density. With the onset of the 
condensate, K~ mesons replace electrons and muons in the core. As a result, K~ mesons 
along with protons maintain the charge neutrality. It was noted that the proton fraction 
became comparable to the neutron fraction in a neutron star including the K~ condensate 
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at higher densities [s|, Q . The appearance of the K ~ condensate would not only influence 
the electron and muon shear viscosities but it will also give rise to a new contribution called 
the proton shear viscosity. 

This paper is organised in the following way. In Sec. II, we describe the calculation of 
shear viscosity in neutron stars involving the K~ condensate. Results are discussed in Sec. 
III. A summary is given in Sec. IV. 

II. FORMALISM 

Here we are interested in calculating the shear viscosity of neutron star matter in the 
presence of an antikaon condensate. We consider neutron star matter undergoing a first 
order phase transition from charge neutral and beta-equilibrated nuclear matter to a K~ 
condensed phase. The nuclear phase is composed of neutrons, protons, electrons and muons 
whereas the antikaon condensed phase is made up of neutrons and protons embedded in the 
Bose-Einstein condensate of K~ mesons along with electrons and muons. Antikaons form 
a s-wave (p = 0) condensation in this case. Therefore, K~ mesons in the condensate do 
not take part in momentum transfer during collisions with other particles. However, the 
condensate influences the proton fraction and equation of state (EOS) which, in turn, might 
have important consequences for the shear viscosity. The starting point for the calculation 



of the shear viscosity is a set of coupled Boltzmann transport equations 
particle species (i= n, p, e, /i) with velocity Vj and distribution function F, 



0,9 



for the ith 



Vi ■ yFi = V Iij. (1) 
The transport equations are coupled through collision integrals given by, 

Ui = (2nhni + 5 l3 ) £ J dPidPi'dPi'W^ , (2) 

where 

T = [F v F f (l - - Fj) - FiFjil - F t ,){l - F f )\ . (3) 

Here p^, Pj are momenta of incident particles and p^, pj> are those of final states. The 
Kronecker delta in Eq. (2) is inserted to avoid double counting for identical particles. Spins 
are denoted by s and Wij is the differential transition rate. The nonequilibrium distribution 
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function for the i-th species F{ is given by 

Ft 

where the equilibrium Fermi-Dirac distribution function /(e^ 



0€i 



(4) 



= e . „. and the de- 

£ i f^i i — i i — i 

1 + e k'F n 

parture from the equilibrium is given by 0. We adopt the following ansatz for 0; [7], [ill] 



1 2 -> -> 

-Ti{viPj - -ViPi8ij){\jiVj + x/jVi - -8ijX7 ■ V) 



(5) 



where tj is the effective relaxation time for the ith species and V is the flow velocity The 
transport equations are linearised and multiplied by (2iih)~ 3 (viPj — ^ViPiSij)dpi. Summing 
over spin Sj and integrating over dpi we obtain a set of relations between effective relaxation 
times and collision frequencies [3] 

{VijTi + P'ifj) = 1 , (6) 



3=n,p,e,fi 

and the effective collision frequencies are 

37r 2 ^ 3 \ f dpidpjdpiidpj 



Vi 



2p 5 F .kTm*( y l + 5, 



(2txK) 



12 



x[^Pi + 3PM'-(Pi-Pi') 2 ] 



(7) 



2*3 



2p F .kTm*Al + 5 ij ) 



E 



dpidpjdpj/dpj 



(2tt^) 



12 



1 „2„2 ^ 2„2 



X [|p 2 Pj' - ^PiPi + (Pi ■ Pj) 2 - (Pi ■ Pj') 2 ] 



The differential transition rate is given by 

^2 w ij = 4(2vr) 4 ^ 2 (5(e i + e,- - - e / )5(p i + - p^ - Pf)Qij , 



(9) 



where Qy =< |tVL,- | 2 > is the squared matrix element summed over final spins and averaged 
over initial spins 

We obtain effective relaxation times for different particle species solving a matrix equation 
that follows from Eq.(6). The matrix equation has the following form: 



V u' np V n j 



\r n / 



(10) 



where, 



v e = v ee + v' ee + u efl + v ep , (11) 
Vy. = + i/ m + + v w , (12) 

V v = Vpp + v' m + V V n + V V e + Vpfi , (13) 
Vn = V nn + V nn + V np . (14) 

It is to be noted here that the proton-proton interaction is made up of contributions from 
electromagnetic and strong interactions. As there is no interference of the electromagnetic 
and strong interaction terms, the differential transition rate for the proton-proton scatter- 
ing is the sum of electromagnetic and strong contributions. This was discussed earlier in 



Ref. 



4j. Therefore, we can write the strong and electromagnetic parts of the effective collision 



frequencies of proton-proton scattering as 



v s +u em 

pp ' pp 5 



/ ts /em /1 r\ 

V VV = Vpp + V pp ■ ( 15 ) 

Here the superscripts 'em' and V denote the electromagnetic and strong interactions. So- 
lutions of Eq. (10) are given below 

Te detA 

Tp detA 

_ {Vp - ~ <jj + Kp ~ ^p)^'p^e - K^'py) + Kp - <pK<e^ ~ U 'pe^) n ^ 

Tn ~ detA [ ' 

where A is the 4x4 matrix of Eq. (10) and detA = [veV^VpVn — v' vn v' nv } — v e i/ 1/ u n — 

^ M ^(^n-^0~^^^" + ^^^ n_I 4 I/ '' I ^e ,/ n]- We obtain r M from r e replacing 
e by fi. In the next paragraphs, we discuss the determination of matrix element squared for 
electromagnetic and strong interactions. 

First we focus on the electromagnetic scattering of charged particles. Here we adopt 
the plasma screening of the interaction due to the exchange of longitudinal and transverse 
plasmons as described in Refs. , [ijjl . The matrix element for the collision of identical 
charged particles is given by M12 = + , where the first term implies the scattering 
channel 12 — > 1'2' and the second term corresponds to that of 12 — > 2'V. The scattering of 
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charged particles in neutron star interiors involves small momentum and energy transfers. 
Consequently both channels contribute equally because the interference term is small in this 
case. The matrix element for nonidentical particles is given by Q, Q] 

47re 2 /jll'j22' jll'.j22' \ 
Mi2-^1'2' = 5— I 2 , -fl2 2 9 / 2 1 TT2 ' (^) 

c 2 \q 2 + Ilf q 2 — cu 2 /c 2 + lLf J 
where q and to are momentum and energy transfers in the neutron star interior. Further 
four-current (Jo, J) and longitudinal and transverse polarization functions (II;, lb) have the 
same expressions as defined in Ref . p| . It is to be noted that J t is the transverse component 
of J with respect to q and the longitudinal component is related to the timelike component 
Jo by the conservation of current |l4j . Polarization functions 11; and lb are associated with 
the plasma screening of charged particles' interactions through the exchange of longitudinal 
and transverse plasmons, respectively. After evaluating the matrix element squared and 
doing the angular and energy integrations, the effective collision frequencies are calculated 
following the prescription of Ref.|7J, |l2j. The collision frequencies of Eqs. (7) and (8) for 
charged particles become 



4 = 4 > ( 18 ) 

where i/]'- and u£- correspond to the charged particle interaction due to the exchange of 
longitudinal and transverse plasmons and is the result of the interference of both in- 
teractions. For small momen ;um and energy transfers, different components of the collision 



frequency are given by [7|, ll2| 



± _ e A a V% (hc\ 1/3 5/3 



e 4 n 2 m*m* 2 



h p%qi 
,im 2e 4 vr 2 m*p|. 



{kTf , (19) 



13 h 4 c 2 p 3 Fi qi 

where i,j = e,fi,p and longitudinal and transverse wave numbers are given by 

2 4 e 2 

h C7T ' J 

?? = £ p% ■ ( 2 °) 

J=e,/J.,p 



The value of a = 2(^) 1 / 3 r(8/3)C(5/3) ~ 6.93 where T(x) and ((x) are gamma and Riemann 
zeta functions, respectively. The shear viscosities of electrons and muons are given by [?| 

riip 2 F Ti 

= -g-r- • (21) 

Here effective masses (m*) of electrons and muons are equal to their corresponding chemical 
potentials because of relativistic effects. It was noted that the shear viscosity was reduced 
due to the inclusion of plasma screening by the exchange of transverse plasmons 
It is worth mentioning here that we extend the calculation of the collision frequencies for 
electrons and muons in Refs. [3, Q] to that of protons due to electromagnetic interaction. 
Before the appearance of the condensate in our calculation, protons may be treated as 
passive scatterers as was earlier done by Ref.p|. However, after the onset of the antikaon 
condensation, electrons and muons are replaced by K~ mesons and proton fraction increases 
rapidly in the system [s|,[lo|. In this situation protons can not be treated as passive scatterers. 

Next we focus on the calculation of collision frequencies of neutron-neutron, proton-proton 
and neutron-proton scatterings due to the strong interaction. The knowledge of nucleon- 
nucleon scattering cross sections are exploited in this calculation. This was first done by 
Ref. 4j. Later recent developments in the calculation of nucleon-nucleon scattering cross 
sections in the Dirac-Brueckner approach were considered for this purpose HQ. Here 
we adopt the same prescription of Ref. [3] for the calculation of collision frequencies due to 
nucleon-nucleon scatterings. The collision frequency for the scattering of identical particles 
under strong interaction is given by 

Vu + v a = - rr? Su , (22) 
3m 2 fi 



iff* 2 /"l fVM 5 ) 10 T 2 /2 

Su = —^- 2 / dx' / dx \ X f2 Q u , (23) 

16ft 7l z Jo Jo v 1 — x 2 — x' 2 

where i = n,p and m n is the bare nucleon mass and Qu is the matrix element squared which 
appears in Eq. (9). Similarly we can write the collision frequency for nonidentical particles 

as 



32m*m* 2 (kT) 2 

3m 2 n h 3 
32mfm*(kT) 2 



_ , „ 

13 - 3m 2 h 3 v ' 



^ = k# ^ • (24) 
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and 



5'. = f- / dx' / cfe± ' v °-^Q ll ■ (25) 



We define x Q = x = x' = tt 2 -, f = V^o (o.25+j () x ) ^ momentum trans- 
fers q = pjr — pj and q' = pj> — p^. We find that the calculation of S^, Su and S'^ 
requires the knowledge of Qu and Qij. The matrix elements squared may be extracted from 
nucleon-nucleon differential cross sections. A detailed discussion on the calculation of matrix 
elements squared from the in-vacuum nucleon-nucleon differential scatterin g cr oss sections 



computed using Dirac-Brueckner approach 



m 



16| can be found in Ref . (l2|, |l3( . We follow 



this procedure in this calculation. It is to be noted here that S pp , S pn and are the new 
results of this calculation. As soon as we know the collision frequencies of nucleon-nucleon 
scatterings due to the strong interaction, we can immediately calculate effective relaxation 
times of neutrons and protons from Eq. (16). This leads to the calculation of the neutron 
and proton shear viscosities as 

2 

n nPF n T n 

^ n 5m* ' 

2 

u pPf p t p , 

* = ■ ( } 

Finally the total shear viscosity is given by 

Vtotal = Vn + Vp + Ve + Vfl • ( 27 ) 

The EOS enters into the calculation of the shear viscosity as an input. We construct 
the EOS within the framework of the relativistic field theoretical model [17I . Here we 
consider a first order phase transition from nuclear matter to K~ condensed matter. We 
adopt the Maxwell construction for the first order phase transition. The constituents of 
matter are neutrons, protons, electrons and muons in both phases and also (anti)kaons 
in the K~ condensed phase. Both phases maintain charge neutrality and (3 equilibrium 
conditions. Baryons and (anti)kaons are interacting with each other and among themselves 



by the exchange of cr, uj and p mesons jjj, 



3 



The baryon-baryon interaction is given by 



S 



the Lagrangian density (l9j, |20] 



B=n,p 



g P Bi^B ■ ft*) 4>B 



+ l - (d,ad»a-mla 2 )-U(a) 



The scalar self-interaction 
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20|,l2l| is 



1 1 

U((t) = - gi m N {goNO-f + - g 2 (g^Ncr) 4 , 



(28) 



(29) 



The effective nucleon mass is given by m* B = m# — g a B<^, where mg is the vacuum baryon 
mass. The Lagrangian density for (anti)kaons in the minimal coupling is given by 0, 1(3, 



22] 



C K = D*KD»K - m%KK 



(30) 



where the covariant derivative is D p = d p + ig^K^^ + W p k^k ' an d the effective mass of 
(anti)kaons is m* K = m K — g a K&- The in-medium energies of K ± mesons are given by 



(31) 



The condensation sets in when the chemical potential of K mesons (fix- — wjt-) is equal 
to the electron chemical potential i.e. \i e = Uk- 



Using the mean field approximation 17|, [18) and solving equations of motion self- 



consistently, we calculate the effective nucleon mass and Fermi momenta of particles at 
different baryon densities. 



III. RESULTS AND DISCUSSIONS 

The knowledge of meson-nucleon and meson-kaon coupling constants are needed for this 
calculation. The nucleon-meson coupling constants determined by reproducing the nuclear 
matter saturation properties such as binding energy E/B = —16.3 MeV, baryon density 
n = 0.153 fm~ 3 , the asymmetry energy coefficient a asy = 32.5 MeV, incompressibility 



K = 300 MeV ,and effective nucleon mass m* N /mN = 0.70, are taken from Ref.|23|]. Next 
we determine the kaon-meson coupling constants using the quark model and isospin counting 
rule. The vector coupling constants are given by 

9wK = T^g^N and g pK = g P N ■ (32) 
9 



The scalar coupling constant is obtained from the real part of K optical potential depth 
at normal nuclear matter density 



Ur (no) 



(33) 



It is known that antikaons experience an attractive potential and kaons have a repulsive 



interaction in nuclear matter 



24 



25 



26 



27 



28 



29J. On the one hand, the analysis of K 



atomic data indicated that the real part of the antikaon optical potential could be as large as 



Ur = —180 ±20 MeV at normal nuclear matter density 24j, |25(]. On the other hand, chirally 



motivated coupled channel models with a self-consistency requirement predicted shallow 



potential depths of -40-60 MeV 
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311 ] . Recently, the double pole structure of A(1405) 



was investigated in connection with the antikaon- nucleon interaction 32j, |33j. Further, the 



highly attractive potential depth of several hundred MeV was obtained in the calculation 



of deeply bound antikaon- nuclear states 3J, |35( . An alternative explanation to the deeply 
bound antikaon-nuclear states was given by others 3(J • This shows that the value of antikaon 
optical potential depth is still a debatable issue. Motivated by the findings of the analysis 
of K~ atomic data, we perform this calculation for an antikaon optical potential depth 
Ur = —160 MeV at normal nuclear matter density. We obtain kaon-scalar meson coupling 
constant g a K = 2.9937 corresponding to UR(n ) = —160 MeV. 

The composition of neutron star matter including the K~ condensate as a function of 
normalised baryon density is shown in Fig. 1. The K~ condensation sets in at 2.43no. Before 
the onset of the condensation, all particle fractions increase with baryon density. In this case, 
the charge neutrality is maintained by protons, electrons and muons. As soon as the antikaon 
condensate is formed, the density of K~ mesons in the condensate rapidly increases and K~ 
mesons replace leptons in the system. The proton density eventually becomes equal to the 
K~ density. The proton density in the presence of the condensate increases significantly 
and may be higher than the neutron density at higher baryon densities [291 ] . This increase 
in the proton fraction in the presence of the K~ condensate might result in an enhancement 
in the proton shear viscosity and appreciable reduction in the electron and muon viscosities 
compared with the case without the condensate. We discuss this in details in the following 
paragraphs. 

Next we focus on the calculation of z/&, and v'-. For the scatterings via the elec- 
tromagnetic interaction, we calculate those quantities using Eqs. (18) and (19). On the 
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other hand, us corresponding to collisions through the strong interaction are estimated us- 
ing Eqs. (22)- (25). In an earlier calculation, the authors considered only S nn and S np |7J] for 
the calculation of the neutron shear viscosity in nucleons-only neutron star matter because 
protons were treated as passive scatterers. It follows from the discussion in the preceding 
paragraph that protons can no longer be treated as passive scatterers because of the large 
proton fraction in the presence of the K~ condensate. Consequently the contributions of 
S pp and S pn have to be taken into account in the calculation of the proton and neutron shear 
viscosities. The expressions of S nn , S pp , S np and S pn given by Eqs. (23) and (25) involve 
matrix elements squared. We note that there is an one to one correspondence between the 

n 

differential cross section and the matrix element squared [13J. We exploit the in- vacuum 

n n 

nucleon-nucleon cross sections of Li and Machleidt [15|, H6J calculated using Bonn interac- 
tion in the Dirac-Brueckner approach for the calculation of matrix elements squared. We 
fit the neutron-proton as well as proton-proton differential cross sections and use them in 
Eqs. (23) and (25) to calculate S nn , S pp , S np and S pn which are functions of neutron (pf„) 
and proton (pf p ) Fermi momenta. The values of pf u ranges from 1.3 to 2.03 fwT 1 whereas 
that of pf p spans the interval 0.35 to 1.73 fm . This corresponds to the density range 
~0.5 to ~ 3.0n . We fit the results of our calculation. Figures 2 and 3 display the variation 
of S nn , S pp , S np and S pn with baryon density. The value of S nn is greater than that of 
S np in the absence of the condensate as evident from Fig. 2. Our results agree well with 
those of Ref. 7|. However, S np rises rapidly with baryon density after the onset of the K~ 
condensation and becomes higher than S nn . It is noted that the effect of the condensate on 
S nn is not significant. Figure 3 shows that S pp drops sharply with increasing baryon density 
and crosses the curve of S pn in the absence of the condensate. However S pp and S pn are 
not influenced by the antikaon condensate. A comparison of Fig. 2 and Fig. 3 reveals that 
Spp is almost one order of magnitude larger than S nn at lower baryon densities. This may 
be attributed to the smaller proton Fermi momentum. We also compute S' np and S' (not 
shown here) and these quantities have negative values. Further we find that the magnitude 
of S' pn is higher than that of S' np . It is to be noted here that S'^ is related to v[j by Eq. 
(24). This is again connected to Eq. (6). Therefore, the values of SL and v'^ can be made 
positive by putting a negative sign between two terms in Eq. (6). 

As soon as we know z/s, we can calculate effective relaxation times using Eq. (16) and 
shear viscosities using Eqs. (21), (26) and (27). First, we discuss the total shear viscosity 
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in nuclear matter without a K~ condensate. This is shown as a function of baryon density 
at a temperature 10 8 K in Fig. 4. Here our results indicated by the solid line are compared 
with the calculation of the total shear viscosity using the EOS of Akmal, Pandharipande 
and Ravenhall (APR) [.'jT/J denoted by the dotted line and also with the results of Flowers 
and Itoh {3, 4|. For the APR case, we exploit the parametrization of the EOS by Heiselberg 
and Hjorth- Jensen [3^]. Further we take density independent nucleon effective masses m* = 
m* = 0.8m n for the calculation with the APR EOS which was earlier discussed by Shternin 
and Yakovlev 112 



]. On the other hand, the results of Flowers and Itoh were parametrized 
by Cutler and Lindblom (CL) [5| and it is shown by the dashed line in Fig. 4. It is evident 
from Fig. 4 that the total shear viscosity in our calculation is significantly higher than other 
cases. This may be attributed to the fact that our EOS is a fully relativistic one. 

We exhibit shear viscosities in the presence of an antikaon condensate as a function of 
baryon density in Fig. 5. This calculation is performed at a temperature 10 8 K. In the 
absence of the K~ condensate, the contribution of the electron shear viscosity to the total 
shear viscosity is the highest. The electron, muon and neutron shear viscosities exceed the 
proton shear viscosity by several orders of magnitude. Further we note that the lepton 
viscosities are greater than the neutron viscosity. On the other hand, we find interesting 
results in the presence of the antikaon condensate. The electron and muon shear viscosities 
decrease very fast after the onset of K~ condensation whereas the proton shear viscosity 
rises in this case. There is almost no change in the neutron shear viscosity. It is interesting 
to note that the proton shear viscosity in the presence of the condensate approaches the 
value of the neutron shear viscosity as baryon density increases. The total shear viscosity 
decreases in the K~ condensed matter due to the sharp drop in the lepton shear viscosities. 
Here the variation of shear viscosities with baryon density is shown up to 3no- The neutron 
and proton shear viscosities in neutron star matter with the K~ condensate might dominate 
over the electron and muon shear viscosities beyond baryon density 3no- Consequently, the 
total shear viscosity would again increase. 

The temperature dependence of the total shear viscosity is shown in Fig. 6. In an earlier 
calculation, electron and muon shear viscosities were determined by collisions only due to 
the exchange of transverse plasmons because this was the dominant contribution [12j . Under 
this approximation, the electron and muon shear viscosities had a temperature dependence 
of T -5 / 3 , whereas, the neutron shear viscosity was proportional to T -2 . The temperature de- 
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pendence of the electron and muon shear viscosities deviated from the standard temperature 
dependence of the shear viscosity of neutron Fermi liquid. However, in this calculation we 
have not made any such approximation. We have considered all the components of effective 
collision frequency which have different temperature dependence as given by Eq. (19). This 
gives rise to a complicated temperature dependence in the calculation of shear viscosity. 
The total shear viscosity is plotted for T = 10 7 , 10 8 , and 10 9 K in Fig. 6. It is noted that 
the shear viscosity increases as temperature decreases. 

The shear viscosity plays an important role in damping the r-mode instability in old and 
accreting neutron stars [39j, |4(J, |4l|, |42|, |43J . The suppression of the instability is achieved by 
the competition of various time scales associated with gravitational radiation {tqr), hyperon 
bulk viscosity (t b ), modified Urea bulk viscosity (t^), and shear viscosity (t S v)- At high 
temperatures the bulk viscosity damp the r-mode instability. As neutron stars cool down, 
the bulk viscosity might not be the dominant damping mechanism. The shear viscosity 
becomes significant in the temperature regime < 10 8 K and might suppress the r-mode 
instability effectively. 

In this calculation, we consider only the antikaon optical potential depth = —160 
MeV. However, this calculation could be performed for other values of antikaon optical 
potential depths. As the magnitude of the K~ potential depth decreases, the threshold of 
the antikaon condensation is shifted to higher densities [9]. On the other hand, hyperons may 
also appear in neutron star matter around 2-3no- Negatively charged hyperons might delay 



the onset of K condensation 



19 



44 



45| . However, it was noted in an earlier calculation 
that £~ hyperons were excluded from the system because of repulsive E-nuclear matter 
interaction and S~ hyperons might appear at very high baryon density [10j. However, the 
appearance of A hyperons could compete with the threshold of K~ condensation. If A 
hyperons appear before K~ condensation, the threshold of K~ condensation is shifted to 
higher baryon density because of softening in the equation of state due to A hyperons. But 
the qualitative results of the shear viscosity discussed above remain the same. 

IV. SUMMARY AND CONCLUSIONS 

We have investigated the shear viscosity in the presence of a K~ condensate. With the 
onset of K~ condensation, electrons and muons are replaced by K~ mesons rapidly. The 
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proton fraction also increases and eventually becomes equal to the neutron fraction in the 
K~ condensed neutron star matter. This has important consequences for the electron, muon 
and proton shear viscosities. We have found that the electron and muon shear viscosities 
drop steeply after the formation of the K~ condensate in neutron stars. On the other hand, 
the proton shear viscosity whose contribution to the total shear viscosity was negligible in 
earlier calculations 0, Q], now becomes significant in the presence of the K~ condensate. The 
proton shear viscosity would exceed the neutron as well as lepton shear viscosities beyond 
3uq. The total viscosity would be dominated by the proton and neutron shear viscosities in 
this case. This calculation may be extended to neutron stars with strong magnetic fields. 

It is worth mentioning here that we adopt the Maxwell construction for the first order 
phase transition in this calculation. Such a construction is justified if the surface tension 
between two phases is quite large [461 ] . Moreover the value of the surface tension between 
the nuclear and antikaon condensed phases or between the hadron and quark phases is not 
iaiown correctly. Therefore, this problem could also be studied using the Gibbs construction 
47]- 

Besides the role of shear viscosity in damping the r-mode instability as well as in pulsar 
glitches and free precession of neutron stars, it has an important contribution in the nucle- 
ation rate of bubbles in first order phase transitions. It was shown earlier that the shear 



viscosity might control the initial growth rate of a bubble 
in connection with antikaon condensation in neutron stars 



49j. This needs further study 
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Fig. 1. Number densities of different particle species as a function of normalised baryon density. 
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Fig. 2. S nn and S np are plotted as a function of normalised baryon density. 
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Fig. 4. Total shear viscosities in nuclear matter without an antikaon condensate corresponding to 
this work (solid line), the parameterization of Cutler and Lindblom (dashed line) and the EOS of 
Akmal, Pandharipande and Ravenhall are shown as a function of normalised baryon density at a 
temperature T = 10 8 K. 
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Fig. 5. The total shear viscosity as well as shear viscosities corresponding to different particle 
species are shown as a function of normalised baryon density at a temperature T = 10 8 K with 
(solid line) and without (dashed line) a K~ condensate. 
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Fig. 6. The total shear viscosity as a function of normalised baryon density at different tempera- 
tures with (solid line) and without (dashed line) a K~ condensate. 
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